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ABSTRACT 


The  one-dimensional,  steady  state  structure  of  a  shock 
vave  in  a  partially  ionized  gas  is  investigated  using  the 
Navier-Stokes  equatiois  for  the  atom,  electron,  and  ion  fluids. 

The  plasma  is  assumed  quasi-neutral  without  a  change  in  ioniza¬ 
tion  across  the  shock,  A  characteristic  feature  common  to 
weakly  and  strongly  ionized  gases  is  the  existence  of  a  broad 
layer  of  elevated  electron  temperature  (thermal  layer)  ahead 
of  the  shock  front  due  to  the  high  electron  thermal  conductivity 
and  of  a  precursor  and  an  imbedded  axial  electric  field  induced 
by  the  charge  separation.  Because  of  the  large  atom-ion 
collision  cross  section  due  to  the  charge  exchange  mechanism, 
the  ion  slip  is  small  when  less  than  30^  of  the  plasma  is  ionized. 
In  a  weakly  ionized  plasma,  the  atom  flow  is  unaffected  by  the 
ionized  particles  and  the  structure  consists  of  an  ordinary  atom 
shock  imbedded  in  the  thermal  layer.  When  the  plasma  is  sub¬ 
stantially  ionized,  the  heavy  particles  are  partially  compressed 
and  heated  in  the  thermal  layer  and  the  ion  and  atom  temperatures 
overshoot  their  downstream  values  in  the  shock.  The  induced 
electric  fields  increase  with  the  degree  of  ionization  and  the 
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free  streaa  Mach  number  vhile  the  shock  thickness  decreases. 

An  experimental  measurement  in  partially  ionized  hydrogen  of 
the  variation  of  potential  in  a  strong  shock  gives  a  shock 
thickness  and  a  potential  rise  across  the  imbedded  shock  in 
good  agreement  vith  the  calculated  values. 

I.  INTRODUCTIOK 

Vhile  the  sttidy  of  a  ihlly  ionized  plasma  is  relevant  to 
sstrophysical  problems,  it  is  more  realistic  to  consider  a 
partially  icmized  plasma  in  problems  of  re-entry  physics  and 
related  laboratory  experiments.  The  purpose  of  this  vork  is  to 
extend  our  solution  for  the  shock  structure  in  a  fully  ionized 
gas^  to  the  case  of  a  partially  ionized  gas,  and  to  investigate 
the  various  plasms  regimes  leading  to  different  types  of  shock  | 
structures.  Ihe  shock  structure  in  an  ionized  gas  is  complex  | 

end  consists  in  general  of  a  thin  region  in  which  most  of  the  ! 

deceleration  and  heating  of  the  flow  occurs,  imbedded  within  a 
much  broader  relaxation  region.  In  our  analysis,  the  term  shock 
refers  to  the  thin  region  while  the  broad  relaxation  region  is 
termed  a  thermal  layer. 

I 

Of  particular  interest  in  the  present  work  is  the  study  of  I 

‘4 

i 

the  effects  of  the  ion-atom  diffusion  or  ion  slip  which  is  measured; 
by  the  velocity  and  temperature  •’.ifferences  between  ions  and  atoms.  ? 
The  negative  electric  field  induced  by  the  charge  separation  in-  ? 
side  the  shock  slows  down  the  ions  and  creates  an  initial  ion  i 


1.  M.  Y.  Jaffrin  and  R.  F.  Probstein,  Phys.  Fluids  7,  Oct.  (19614). 
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slip,  which  vanishes  at  the  end  of  the  shock  after  a  sufficient 
number  of  ion-atom  collisions. 

As  in  our  previous  work,  we  seek  the  solution  for  a  steady 
one-dimensional  shock  with  no  applied  external  electric  or 
magnetic  fields.  The  plasma  is  assumed  here  to  be  a  ternary 
mixtxxre  of  rtoms,  singly  ionised  ions,  and  electrons.  The 
analysis  is  based  on  a  three-fluid  continuum  approach  using  the 
Uavier-Stokes  equations  as  a  model. 

In  order  to  make  the  problem  mathematically  tractable,  we 
introduce  three  basic  assumptions: 

a.  Frozen  ionization;  The  increase  in  temperature 
across  the  shock  produces  further  ionization  of  the  gas  following 
the  shock,  however,  it  is  assumed  that  the  characteristic  chemical 
lengths  are  large  compared  to  the  shock  thickness  so  that  no 
ionization  or  recombination  reactions  take  place  in  the  shock 
itself, 

b.  Quasi-charge  neutrality;  It  is  assumed  that  the 
over  all  charge  separation  is  negligible,  a  condition  satisfied 
when  the  Debye  length  is  much  smaller  than  the  mean  free  path, 
which  is  met  for  most  plasmas. 

c.  Small  ion  slip;  This  follows  from  the  fact  that 

because  of  the  large  ion-atom  cross  section  due  to  the  charge 

exchange  mechanism,  the  ion-atom  mean  free  path  i.  ,1s  in  general 

2 

small  compared  to  the  shock  thickness  As  a  result  the  ion 

slip,  which  is  of  the  order  of  t.  /A  ,  is  small  and  me^  be 

1^2  s 

computed  by  an  Iterative  method. 


The  shock  structure  in  a  fully  ionized  gas  has  been  in- 

1-li 

vestigated  by  several  authors  but  the  case  of  a  partially 
ionized  gas  has  received  less  attention.  Grewal  and  Talbot^, 
assuming  a  Mott-Smith  solution  for  the  ions  and  atoms,  and 
neglecting  the  ion  slip  found  the  existence  of  a  broad  thermal 
layer  in  front  of  the  shock  where  the  electron  temperature 
rises.  But  they  did  not  investigate  the  structure  of  the  ion- 
atom  shock  and  the  assumption  of  small  electron  energy  restricts 
their  solution  to  the  weakly  ionized  case,  Pikel'ner^  con¬ 
sidered  the  inviscid  transition  region  between  a  magnetohydro¬ 
dynamic  shock  and  the  downstream  state  but  he  did  not  take  into 
account  the  charge  exchange  mechanism  which,  by  considerably 
increasing  the  atom-ion  cross  section,  strongly  inhibits  the 
ion  slip  in  the  relaxation  region. 

II.  BASIC  EQUATIONS 

Ibere  are  for  each  type  of  particle,  atom,  ion,  and  electron 
the  three  continuum  equations  of  mass,  momentum,  and  energy.  The 
electric  field  whose  only  non-zero  component  is  is  governed 
by  Poisson’s  equation 


2.  J.  D,  Jukes,  J.  Fluid  Mech.  3,  275  (1957) 

3.  V.  D,  Shafranov,  Zh.  Eksperim,  i  Teor,  Fiz.  32,  1^53  (1957) 

(Eng.  Transit  Soviet  Phys.  JETP  5,  1183  (1957)] 

1|.  0,  W.  Greenberg,  H.  K,  Sen,  and  Y,  M.  Treve,  Phys.  Fluids  3, 

379  (I960) 

5.  M.  S.  Grewal  and  L.  Talbot,  J.  Fluid  Mech.  3^,  573  (1963) 

6,  S.  B.  Pikel’ner,  Zh,  Eksperim,  i  Teor.  Fiz.  36,  1536  (1959) 

[Eng.  Transl.t  Soviet  Phys.  JETP  7,  9  (195971 
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dE' 

=  l*Tte(nI  -  n' ) 
dx'  1  e 


(2.1) 


Here  e  is  ,the  raafmitude  of  the  electronic  charge,  x'  the  axial 
distance,  and  n'  the  number  density;  the  subscripts  i,  e  and 
a,  denote  ions,  electrons  and  atoms  respectively  and  the 
"primes"  are  used  to  distinguish  the  physical  variables  from 
the  dimensionless  variables  which  will  be  used  later. 

The  net  current  vanishes  so  that  we  have 


J  =  0  =  (n!u!  -  n 'u* )  e 
''x  1  i  e  e 


or 


where  u’  is  the  velocity.  On  integrating  the  steady  one  dimensional 
mass  conservation  equations  for  the  ions  amd  electrons,  when  there 
is  no  ionization  or  recombination,  the  constants  of  integration 
and  are  the  same,  giving: 

n!u!  =  n'u*  =  C.  =  C  .  (2.2) 

i  i  e  e  i  e 

The  conservation  of  mass  for  atoms  yields 


n’u’  =  C 
a  a  £ 


(2.3) 
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The  constants  of  integration  C.  and  C  are  related  to  each  other 

& 

through  the  degree  of  ionization  o.  In  the  undistiu’bed  plasisa, 
all  three  species  have  the  sene  nean  velocity  and  the  degree  of 
ionization  a  can  be  expressed  by  the  relation: 


n!  +  n’  C.  +  C 
1  a  i  a 


The  conservation  of  noaentua  for  the  three  species  yields 
the  following  equations: 


du*  dp*  du* 

j  j  dx*  dx*  dx*  dx*  J  J  x  Jk 


(2.5) 


where  m  is  the  particle  mass,  y"  the  longitudinal  coefficient 
of  viscosity,  p*  the  pressure  and  P..  represents  the  longitudinal 

jK 

momentuH  transfer  between  species  J  and  k  and  where  J  =  i,  e,  a 
with  the  convention  e.  «  +e,  e  *  -e,  e  =0,  The  summation 
is  made  for  k#  J,  The  longitudinal  coefficient  of  viscosity  y* 
is  related  to  the  claissical  coefficient  of  viscosity  y  by  the 
relation 


y”  =  l«/3  y  . 


The  corresponding  energy  equations  are 


du m , 


5  d  (p.^. )  .  ^  -4)  -  A  (u*y”  -4)  n'u* 

2  ^  J  dx^  ^  ”  2  "j  J  dx' 


=  I  (CjK  .  ujPj,)  , 
J 


(2.6) 
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where  k  is  the  thernal  conductivity,  T'the  absolute  tenperatare 
and  denotes  the  kinetic  energy  transfer  froa  the  k- 
particles  to  the  J-particles.  The  basic  systen  is  cotapleted 
by  the  equations  of  state  for  the  three  species 


pj  =  n-  kT-  .  (2.7) 

With  the  dissipative  coefficients  and  transfer  terns  assmaed 

to  be  known  functions  of  the  dependent  variables,  then  Eqs. 

(2.1)  to  (2.7)  constitute  a  system  of  thirteen  equations  for 

the  thirteen  mknowns  u',  n*,  T'  p*  and  E*.  However,  since  the 

J  J  o  ^ 

continuity  and  state  equations  are  algebraic,  and  six  in  number, 
the  problem  is  reducible  to  the  solution  of  seven  simultaneous 
differential  equations  for  the  appropriate  seven  unknowns. 

We  may  simplify  the  system  of  equations  first  by  noting 
that  from  Hewton*s  third  law  of  motion. 


H  0 


(2.8) 


4  U’.P 


J  JK 


(2.9) 


Secondly,  we  note  that  the  electron  inertia  and  viscous  terms  are 
much  smaller  than  the  corresponding  ion  and  atom  terms  and  may 


therefore  be  neglected. 
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After  oqployinc  the  aboTe  conditions  and  approxisaticms  the 
basic  set  of  shock  structure  equations  vill  consist  of  the  tvo 
algebraic  continuity  equations  (2.2)  and  (2.3),  Poisson's 
eqxiation  (2.1),  and  the  folloving  six  differential  equations: 

a.  The  electrtxi  aosjentun  equation  (2,5),  J  =  e 


(2.12) 
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Eere  the  atoa  nass  is  assised  to  be  equal  to  tne  ion  nass  and 
the  constant  of  integration  ?  is  detemined  by  the  boundary 
ccHiditions. 

d.  Ihe  electron  internal  energy  equation,  which  is 
obtained  by  nuitiplying  the  electron  conentun  equation  by  u^ 
and  subtracting  the  result  froa  the  electron  energy  equation. 
It  Bay  be  written 


,  dT'  kT»  du*  ,  dT* 

—  C  k  — +  C  ^  ^  ^  ( Y- 

2  1  dx’  i  u*  dx’  dx*  ^  e  dx*^ 

e 


= 


ea 


^ei 


(2.13) 


e.  The  ion  internal  energy  equation,  obtained  by  the 
same  process  as  Eq.  (2.13),  which  may  be  written 


dT! 


kT!  du! 


"“1.2  d  ,  "n, 


J  i  ^  _ i  _ i  tf  / _ Q  /  _ i% 

2  1  dx’  1  u!  dx'  1  dx'  dx'  i  dx' 


=  +  C. 

la  xe 


(2.1J4) 


f.  Ihe  plasma  energy  equation,  obtained  by  adding 
Eqs.  (2.6)  for  the  three  species  with  the  aid  of  Eq,  (2.9)  and 
integrating  the  resulting  equation  to  give 


C. 


dT' 
a 


m.  ~  u’^  +  m,  u!^  +  I-  klC  T’  +  C.{T!  +  T’ ) )  -  ic  -r~7 
i2  a  i2  i  2  aa  l  i  e'  a  tlx* 


dT!  dT'  du'  du! 

i  ®  #  ft  ^  t  ^  tl 

-  tc.  -  IC  -5-^  -  u*u  -T— r  •  u:u,  -r-r  ="  H 
1  dx*  e  dx'  a  a  ax'  i  i  dx' 


9 


(2.15) 


the  constant  of  integreticn  H  is  evaluated  at  the  boundary 


conditions. 

III.  DISSIPATION  AUD  TRAHSFER  COEFFICIENTS 
We  nov  evaluate  the  dissipaticm  and  transfer  coefficients 
which  are  present  in  the  shock  structure  equations  (2.1)  -  (2.15). 
'3ie  classical  coefficients  of  viscosity  and  thensal  conduction 

7 

for  a  pure  nonatoBiic  gas  are  given  respectively  by  the  relations 
u  *  n*iBc(^)  ,  (3.1) 

1/2 

where  c  «  (8  kT*/ina)  is  the  mean  thermal  speed,  v  the  col¬ 
lision  frequency  and  k  the  Boltzmann  constant.  The  collision 
frequency  for  a  p’ore  gas  is  equal  to: 

V  =  »^  n*Q'  c  ,  (3.3) 

aa  ’ 

where  Q'  is  the  effective  hard  sphere  cross  section  for  atom- 
atom  collisions. 

The  relations  (3.1)  and  (3.2)  can  be  extended  to  the  con- 

Q 

ponent  I  of  a  gw  mixture  as  follows  : 


7.  S.  Chapman  and  T.  G.  Cowling  The  Mathematical  TheoiT  of  Non- 
Uniform  Gases,  (Cambridge  lAiiversity  Press,  Hew  York,  1952) 
pp,  lOO-lQl), 

8.  J.  A.  Fay,  "Hypersonic  Heat  Transfer  in  the  Air  I^aninar  Boundary 
Layer,"  Avco  Everett  Research  Lab.  AMP  71,  March  1962. 
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es  M 


<.  - 


15k 


(3.^»b) 


where  is  the  total  collision  frequency  for  momentuEi  exchange 
given  by 


=  I 

k 


(3.5) 


Here,  k  is  sunmed  over  all  components  of  the  mixture, 

2  2  1/2 

m^j^  =  ^  ®k^  reduced  mass,  ^ 

mean  relative  thermal  velocity  and  is  the  fraction  of 

momentum  transferred  per  collision  between  2  particles  of 
arbitrary  mass.  The  quantity  QJ  is  the  collision  cross  section 
for  particles  i  and  k.  For  a  pure  gas,  reduces  to  v  defined 
by  Eq.  (3.3). 

On  developing  the  relation  (3.5) t  the  atom,  ion,  and  electron 
collision  frequencies  may  be  written  respectively  as 


kT*  ,  •  T»  +  Tf  ,  8kT*  -  2m 

'’a  =  ^  n^o;  [8k(  n’O*  (-~)^^^ 

a  a  aa  Tim.  1  ia  Tim.  e  e*  nm  m. 

i  i  e  i 


(3.6) 


T'  +  T!  .  kT!  ,  kT'  ,  2m 

V  =  n'Q!  tek(  “  *  nJQ!  (16  *  n-Q’  (8  -i 

1  a  la  «m.  i  ii  Tim,  e  el  Tira  m, 

i  i  e  i 


(3.7) 
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V 

e 


^  lo  ^  lo  ^  to 

2n'Q'  (8  ^  2ii!QV(8  +  n'Q»  (l6  — 

a  ea  to  i  ei  to  e  ee  to 

e  e  e 


C3.8) 


’Hie  last  terns  of  Eqs.  (3.6)  and  (3.7)  represent  the  effect  of 
the  electron-heavy  particle  collisions  and  are  of  the  oi.'er  of 
where 


c 


n. 

1 


1 

275 


for  Argon. 


These  terns  are  very  snail  and  will  be  neglected  in  comparison 
to  one.  It  is  now  necessary  to  evaluate  the  various  cross 
sections  remaining  in  the  above  relations. 

The  cross  sections  for  collisions  between  electrically 
charged  particles  are  deduced  from  the  value  of  the  viscosity 

9 

coefficient  for  an  ionized  gas  given  by  Chapman  and  Cowling  . 
Using  the  Spitzer^^  result  that  the  Debye  length  is  the  proper 
cutoff  impact  parameter,  one  finds 


U 


>  5  /ina/2 

-  E 


e  logg  A 


(3.9) 


where  A  is  a  dimensionless  cutoff  impact  parameter  equal  to 
3(k^T'^/im’ )^*^^/2e^.  Comparing  Eqs.  (3.1)  and  (3.9)  we  obtain 
the  following  values  for  the  Coulomb  cross  section; 

ne^  log  A 

Q:  =  - ,  J  =  i,  e  .  (3.10) 

2(kTp 

9.  See  Ref.  7,  p.  179. 

10.  L.  Spitzer  Jr. ,  Physics  of  Fully  Ionized  Gases  ( Interscience 
Publishers,  Inc.,  Hew  York,  195d)  p.  72. 
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Since  electrcm  temperature  is  the  relevant 

temperature  ,in  the  calculation  of  ion-electron  collision  cross 
section  and 


QV  =  Q* 

ei  ee 


(3.11) 


The  values  of  the  other  required  cross  sections  Q*  ,  Q!  , 

and  Q’  depend  on  the  gas  considered.  They  will  here  be 

evaluated  from  experimental  results  for  Argon.  Ihe  atom-atom 

collision  cross  section  Q*  is  obtained  from  the  values  of  the 

aa 

viscosity  coefficient  y  given  by  Amdur  and  Mason^^.  At  high 

Cl 

-7  3  A 

temperatures  u  'v-  31.10  T*  g/cm  sec,  which  corresponds  to 
a  ^  a 


Q'  =  170.10“^^/T*^^^  cm^  . 
aa  a 


g 

Experimental  data  ccmipiled  by  Fay  show  that  the  atom-ion  cross 
section  Q*^  is  much  bigger  than  the  atom-atom  cross  section  be¬ 
cause  of  the  charge  exchange  mechanism.  This  cross  section 
decreases  very  slowly  with  the  temperature  and  will  be  taken 
constant  and  given  by 


1  p 

Q!  =  ll^O  •  10  cm 
la 


11,  I.  Amdur  and  E.  A.  Mason,  Phys.  Fluids  1^,  370,  (1958) 


12  13 

Ihe  atoei--clcctron  cross  section  Q*  *  has  a  lov  nininim  of 

ea 

0.25  •  lO”^^  cx?  at  =  J*600°K  because  of  the  strong  Ransauer 
effect.  It  increases  rapidly  for  higher  temperatures  up  to  10^  °K 
and  decreases  again.  For  a  range  of  temperature  up  to  5  *  10^  °K 
we  shall  use  the  following  approximate  relations  obtained  from  Ref.  12  by 
curve  fitting 


Q’  »  (-0.35  +  0.775  lO"**  T*)  x  lo"^^  cm^  T»  >  10^  °K 
ea  e  e 

Q*  =  (0.39  -  0.551  10"^  T*  +  0.595  lO"®  cn^  T*  <  10 

ea  e  e  e 


In  any  case  the  momentum  and  energy  exchange  between  electrons  and 

atoms  is  so  slow  that  the  structure  of  the  shock  itself  depends 

only  very  weakly  on  the  cross  section  Q'  , 

68. 

For  generality  the  analysis  will  be  carried  out  with  the  cross 

sections  expressed  as  known  functions  of  the  temperature  which  are 

not  explicitly  si)ecified.  The  results  obtained  will  be  qualitatively 

valid  for  other  values  of  the  cross  sections  provided  that  the 

essential  feature  Q!  >>  Q*  is  retained, 

ia  aa 

The  viscosity  and  thermal  conduction  coefficients  for  the  atom 
gas  may  then  be  written  in  the  form 


u 


a 


(mn.  kT* ) 
1  a 


1/2 


[1  + 


n]  q; 


n'  Q’ 

a  aa 


rpf  x  ipi 

ia  ,_a _ i.\l/2') 

'  OfPI  f  J 


2T« 

a 


(3.12a) 


12.  I.  P.  Shkarofsky,  M,  P.  Bachynski  and  T.  W,  Johnston,  Planetary 
and  Space  Science  2^*,  1961 

13.  S,  C.  Brown,  Basic  Data  of  Plasma  Physics  (The  Technology  Press  of 
the  Massachusetts  Institute  of  Technology,  and  J.  Wiley  and  Sons 
Inc.,  N,  Y.,  1959)  p.  19. 
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ti 


Q* 

^ia 


T* 

r  a 


oa 


2^ 

a 


(3.12b) 


We  note  that  the  viscosity  and  thenaal  conduction  coefficients 
of  the  atom  gas  are  reduced  by  the  presence  of  the  ions. 

For  the  ion  gas  these  coefficients  are  given  by 


5n! 


iSq 


n 


r  T! 


ia  a 


T*  +  T! 


1)1/2 


11 


n"^ 
a  ia 


2T! 


TJ 


,)l/2] 


(3.13b) 


K.  = 


wk  1/2 

Aq’  n*  m,  i 

ia  a  i 


*  '^1  1/2  ”i  ‘^ii 

'-V^>  '  < 

a  ia 


F 

a 


j1/23 

nTF 


(3.13b) 


while  for  the  electron  gas 


K 

e 


1/2 

T>i~)  '  k 

PI 

e 


/2  n'  Q' 

6l4Q»  (1  +  »^)  (1  +  ® 


ee 


(1  +  ^)  n'  Q' 
e  ee 


(3.1M 


“l6“ 


Ihe  Yiscosity  and  ther»al  conduction  coefficients  of  the  ionized 

particles  are  here  reduced  by  the  presence  of  atoms ,  The  electron 

viscosity  coefficient  is  smaller  than  by  a  factor  of  the  order 
®e  1/2 

of  £  »  (— )  and  the  electron  viscosity  term  may  be  neglected 

®i 

in  the  structure  equations,  as  it  was  noted  earlier. 

From  Eqs,  (2,12)  and  (2.15),  we  note  that  for  given  values  of 
the  Mach  number,  velocities  and  temperatures,  the  velocity  and 
temperature  gradients  are  inversely  proportional  to  the  viscosity 
and  the  thermal  conduction  coefficients.  Since,  from  Eq,  (2,10) 
the  electric  field  is  proportional  to  a  temperature  and  velocity 
gradient,  the  whole  shock  structure  depends  critically  on  the 
behavior  of  the  dissipation  coefficients.  As  the  viscosity  and 
thermal  conduction  of  the  atoms  coefficients  are  progressively  re¬ 
placed  by  those  of  the  ions  as  the  main  dissipative  mechanisms, 
the  plasma  regime  changes  from  weakly  ionized  to  fully  ionized  and 
the  shock  structure  varies  accordingly,  Ihe  different  plasma  regimes 
for  Argon  are  illustrated  in  the  o,  T*  diagram  of  Fig.  1,  When 
n!Q!.  <<  n’Q!  and  njQ*  «  n’Q'  ,  it  may  be  seen  from  Eqs,  (3,12) 

i  xx  &  xft  X  xft  8t  ftCl 

and  (3,13)  that  the  ion  and  atom  dissipative  coefficients  depend 
mostly  on  the  atoms.  Moreover,  if  a  <<  1  the  ion  dissipation  is 
negligible  compared  with  the  atom  dissipation  and  the  region  delimited 
by  the  three  preceding  inequalities  represents  the  weakly  ionized 
regime.  On  the  other  hand,  when  >>  and  "a^la* 

the  atom  collisional  dissipation  is  much  smaller  than  the  ion 
dissipation  and  the  ion  viscosity  and  thermal  conduction  approach 
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tbeir  fully  ionized  value:  this  we  term  the  quasi-fully  ionized 

_3 

region.  Finally,  the  intermediate  range  between  a  =  10  and 
o  =»  0,5  represents  a  mixed  regime  that  we  denote  by  partially 
ionized. 

The  momentum  exchanges  between  the  different  species  of  a 

ih 

partially  ionized  gas  have  been  computed  by  Zhdanov  ,  On 
neglecting  thermal  diffusion  but  extending  Zhdanov's  results  to 
a  multi-temperature  model,  the  momentum  transfer  per  imit  volume 
from  species  k  to  species  j  is  given  by 


where 


I,  8kT;  8kr  , 

^Jk  "  3  Sk  • 


(3.16) 


The  quantity  Z  is  of  the  same  order  as  the  collision  time.  In 
jK 

particular,  the  ion-atom  momentum  transfer  is 


'■ai  =  -  I  -  “i’  •  <3.17) 


T'  T! 

6  i 

Using  the  inequality  —  »  —  the  electron-atom  momentum  transfer  is, 

m  m. 

e  i 


Q  m  kT'  - 

P  =  -  i  n'u'{-^ - Q*  (u'  -  u') 

ae  3  aew  eaa  e 


(3.18) 


ll4.  Zhdanov  V,  M. ,  PMM  26 ,  280,  (1962),  [Eng.  Transl.:  Applied 
Mathematics  and  Mechanics  ^iOl  (1962).] 
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Since 


1/2  1/2 


a.  »  K 


it  follovs  that 
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Ihe  ic»-electron  BOBentios  transfer  is  given  by 


m  kT» 


P^l  .  -  I  ^  .  (3.19) 


The  total  ion-electron  energy  transfer  has  been  conputed  by 

Burgers"^^,  Expanding  Burgers’  result  in  powers  of  the  small 
°e  1/2 

quantity  e  »  (— -)  we  find  that  the  ion-electron  energy 

1 

transfer  may  be  written 
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(3.20) 


Petschek  and  Byron^^  have  computed  the  atom-electron  energy 
transfer  due  to  the  random  motion  of  the  particles  for  arbitrary 
interaction  laws  which  is  given  by 


15,  J,  M.  Burgers,  in  Plasma  I^amics,  edited  by  F.  H.  Clauser 
(Addison-Wealey  Publishrng  Co.^  i^eading.  Mass.,  I960)  p.  156. 

16.  H.  Petschek  and  S.  Byron,  Annals  of  Physics  1^,  270,  (1957). 
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where  c  is  the  randoa  electron  velocity,  a(c  )  the  cross 
“  e 

section  and  f  the  Maxwell  distribution  function.  The  in- 
e 

tegral  tern  in  Eq,  (3.21)  represents  the  nean  value  of  the 

3 

product  c  o(c  ).  Recalling  that  Q*  is  the  cross  section 
e  e  ^ea 

averaged  over  a  Maxwell  distribution  and  approxinating  the 
nean  value  of  the  product  c^  product  of  the 

averages  of  the  two  terras,  we  obtain 
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The  total  energy  transfer  is  the  sun  of  the  energy  transfers 
due  to  the  random  and  directed  motions  of  the  particles  and  is 
given  by 
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The  ion -atom  energy  transfer  is  derived  in  Appendix  A  and 
the  result  is 
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IV. 

At  upstreiia  infinity  (x  -*■  -  •)  and  at  dovnstreaa  infinity 
(x  -♦  +  •)  all  gradients  vanish  and  the  systeo  of  equations  (2.1)  - 
(2.15)  reduces  to  the  Hugoniot  conditions.  Ve  note  that  on  each 
side  of  the  shock,  the  three  species  have  the  sane  velocity  and 
teaperature,  and  the  electric  field  vanishes. 

Since  the  shock  thickness  depends  generally  on  the  dovnstrean 
state,  we  choose  to  nomalize  the  velocities  and  teaperatures  with 
respect  to  their  dovnstrean  values  after  the  shock.  In  what  follows 
VC  denote  the  upstrean  state  by  subscript  1  and  the  downstrean 
state  by  subscript  2.  On  this  basis  we  set 
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vhere  is  the  actual  potential. 

We  also  introduce  the  diaensionless  charge  separation  defined 
the  relation 


n!  -  n* 
1  e 


i^.U) 


lie  independent  variable  is  cade  diaensionless  with  respect  to 

the  shock  thickness  A  by  the  relation 

s  ^ 


X  =  x'A  , 
s  s 


{h.5) 


Ihe  cross  sections  Q*  are  noraalized  with  respect  to  their 

ao 

downstreaa  values  by  the  relation 


Q  =  Q'  /Q' 

a6  aB, 


(1*.6) 


It  is  convenient  at  this  .point  to  introduce  the  following 
mean  free  paths  which  are  related  to  the  various  downstream  cross- 
section  Q'  „  : 

atom-atom  mean  free  path 


t  » 


"I  Q’  . 


a2  aa^ 


(»4.7a) 


ion-ion  mean  free  path 
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atoB-ion  oean  free  path 
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GtMi^clectron  mean  free  path 
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Since  n!  =  n*  ,  Q*  “  Q*i  *  QI .  ,  the  ion-electron  and  the 

Ig  .2  “2  '^2  “2 

electron-electron  mean  free  paths  are  both  equal  to  the  ion- 
ion  mean  free  path. 

Finally,  we  define  the  downstream  plasma  I-Iach  number  by 
the  relation 
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where  as  previously  the  electron  mass  is  neglected  in  cwnparison 
to  the  ion  mass. 

If  we  emplov  the  non-dimensionalization  of  Eqs.  (^<.1)  -  {^<.8) 
neglecting  terms  of  order  e  and  smaller,  such  as  the  electron 
viscosity  and  inertia  terms,  then  we  may  write  in  dimensionless  form: 

Poisson's  equation  (from  Eq,  (2.1)) 
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On  the  other  hand,  the  electron  nonentun  equation  is  given  fraa 
Eq.  (2.5)  with  J  =  e  by 
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From  the  above  relation  we  see  that 

^“2 

E  =  0(— )  , 

s 

while  from  Eq.  (^.9) 
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In  most  plasmas,  the  Debye  length  is  much  smaller  than  the  mean 
free  path,  and  therefore  much  smaller  than  the  shock  thickness. 

It  follows  that  the  dimensionless  quantities  E  and  6  are  very  small 
and  can  be  neglected  in  comparison  to  one.  Therefore,  in  all 
structure  equations,  except  Poisson's  equation  {h,9)  we  set 
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We  note  that  the  ion  and  electron  fluids  nove  together  although 
their  teeperatures  nay  be  different. 

With  the  approxiaiaticHis  noted  above  our  basic  dinensionless 
equations  in  addition  to  Poisson's  equation  nay  be  written: 
Electron  monentun  equation  (froo  Eq,  (2.5)) 
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Ion-electron  Homentun  equatir-.  (from  Eq,  (2.11)) 
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Plasna  mooentun  equation  (from  Eq.  (2.12)) 
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Electron  internal  energy  equation  (from  Eq.  (2.13)) 
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Ion  internal  energy  equation  (froa  Eq,  (2.lU)) 
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Plasaa  energy  equation  (froa  Eq.  (2,15)) 
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Bie  diHCDsionless  potential  nay  be  obtained  from  the  relation 

2 

The  quantities,  a,  b,  c,  d,  f,  F  ,  F. ,  and  F  appearing  in 
the  dimensionless  equations  (1»,10)  -  (1*,15)  are  defined  as  follows: 
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When  the  quantities  F^,  F^,  and  F^  are  replaced  by  the 
definitions  given  above,  Eqs.  (^,9)  -  (^tl5)  represent  a  set  of 
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seven  equations  for  the  seven  unknown  variables  u  ,  u. ,  T  ,  , 

a*  i*  a*  i* 

T  ,  E  and  6  eis  a  function  of  the  independent  variable  x  .  Had 

G  S 

we  normalized  the  variables  with  respect  to  their  upstream  vsilues 
we  would  havp  received  precisely  the  sane  set  of  equations  with 
the  subscript  1  replacing  the  subscript  2  in  the  parameters. 

Any  state  of  the  plasma  will  be  represented  by  a  point  in 
the  seven  dimensional  phase  space  whose  coordinates  are  the 
dependent  variables.  The  states  1  and  2  upstream  and  downstream 
of  the  shock  respectively,  are  singular  points  of  the  system  of 
Eqs.  (^.9)  -  (*♦.15).  The  dimensionless  Hugoniot  relations  obtained 
by  setting  the  derivatives  equal  to  zero  specify  the  coordinates 
of  these  singular  points.  In  particular  the  coordinates  of  point 
1  are 
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while  the  coordinates  of  point  2  are 
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Here  the  Mach  nunber  behind  the  shock  is  related  to  the  free 
strean  Mach  number  by  the  relation 


(^♦.23) 


It  is  important  to  note  that  the  Hugoniot  conditions  depend  on 
the  plasma  Mach  nxanber  only.  For  the  same  velocity  and  tem¬ 
perature  conditions  upstream,  the  plasma  Itech  number  and  therefore 
the  shock  strength  are  reduced  when  the  degree  of  ionization  in¬ 
creases, 

A  shock  will  exist  whenever  u^  >  1,  which  requires  >  1, 

The  atom  and  ion  Mach  numbers  are  greater  than  the  plasma  Mach 

1/2 

number  by  the  factor  (l  +  a)  ,  and  the  atom  and  ion  flows  are 

supersonic  when  >  1.  The  electron  Mach  number  is  proportional 

to  e  times  the  plasma  Mach  number,  so  that  the  electron  flow  is 

subsonic  unless  M^  is  very  large. 

The  shock  thicknesses  A  in  which  the  various  collisional 

s 

effects  are  important  can  now  be  estimated  from  the  appropriate 

terms  in  the  dimensionless  equations.  For  example,  we  note  from 

Eq.  (^,13)  that  the  electron  thermal  conduction  term  is  of  the 

order  of  £,.  /eA  ;  it  is  therefore  of  order  one  when  A  is 
ii^  s’  s 

0(£.  /e).  This  implies  that  the  electron  thermal  conduction 

can  broaden  the  shock  to  a  thickness  of  the  order  of  £.,  /c.  In 
the  weakly  and  partially  ionized  regimes  where  the  dimensionless 
expressions  F  ,  F. ,  and  F  are  0(l),  the  possible  shock  thicknesses 
and  their  corresponding  collisional  mechanisms  are: 
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Mechanisms 

Electron  thermal  conduction. 


Shock  thickness 


I  /e 
ea^ 


aa^ 

< 

1. 

la^ 

c 

£. 

la^ 


electron-ion  energy  transfer, 

electron-atom  energy  and  momentimi  transfer, 

atom  thermal  conduction  and  viscosity, 

ion  thermal  conduction  and  viscosity, 

atom-ion  momenttun  and  energy  transfer. 

Obviously  the  shock  structure  will  depend  on  the  relative 

magnitude  of  the  various  mean  free  paths.  The  values  of  the 

ratios  £.  /£  ,  £  /£..  and  £.  /£..  are  plotted  in  Fig,  2 

xa^  aa^  ^^2  ^^2  ^^2  ^^2 

as  a  function  of  the  temperature  for  different  values  of  the 

degree  of  ionization  for  Argon.  It  may  be  noted  that  the  ratio 

£.  /£  is  independent  of  a  and  is  always  small  compared  to  one. 

1  ^^2 

The  ion-ion  mean  free  path  £  .  is  larger  than  the  atom-atom  mean 

free  path  £  for  low  degrees  of  ionization  but  it  is  smaller  than 
aa2 

£  for  a  =  0.5  and  T'<  2.5  10^  °K. 
aa^ 


In  the  quasi-fully  ionized  regime  F  =  0{£  /£.  ), 

s  ^ 

F.  »  0(£,  /£..  )  and  the  atom  thermal  conductivity  and  viscosity 

lead  to  a  shock  thickness  0(£.  )  while  the  ion  thermal  conduction 

and  viscosity  correspond  to  a  shock  thickness  0{£..  ),  The  shock 

1I2 

thidmesaes  corresponding  to  the  other  collisioned  mechanisms  are  not 
modified.  However  we  shall  not  consider  this  regime  because  of  its 
similarity  with  the  fully  ionized  case^. 
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V.  WEAKLY  IONIZED  PLASMA  g  «  1 

We  consider  first  the  case  of  the  weakly  ionized  plasma, 

because  of  the  simplifications  introduced  by  the  assumption 

a  «  1,  Ihe  atom  gas  is  indeed  unaffected  by  the  charged 

particles  provided  that  the  degree  of  ionization  is  very  small. 

We  note  that  the  terms  representing  the  charged  particles  in  Eqs. 

(h,12)  and  (1*,15)  are  proportional  to  a  and  can  be  neglected  when 

a  <<  1,  Moreover,  when  a  is  so  small  that  aZ  1 1,  «  1,  for 

aa^  ia^ 

example  when  a  -  lO”  ,  Eq.  (ii.l8)  yields 


F  'v-  Q 
a  aa 


(5.1) 


With  the  approximations  described  above,  Eqs,  (^.12)  and 
(^.15)  govern  the  atom  velocity  and  temperature  distributions 
in  the  atom  shock  and  are  completely  uncoupled  from  the  other 
variables. 

We  first  examine  the  thermal  layer  of  thickness  /c. 
Should  a  discontinuous  solution  arise  within  this  layer,  we  in¬ 
troduce  at  the  discontinuity  a  shock  layer  whose  thickness  is  the 
next  smaller  length  until  a  continuous  transition  is  found  from 
the  upstream  state  to  the  downstream  one. 

When  A  =£..  /c,x  is  replaced  by  x*  =  cx'/£..  ,  On 
8  ii^  *8  ^  •'  ii^ 

neglecting  terms  of  order  of  a  and  c  in  Eqs,  (^,9)  -  (^.l6),  we 
obtain  a  simplified  system  of  shock  layer  equations.  From  Eqs, 
{U,12)  and  (1|.15) 
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T 

.  1)  +  1  =  0  ,  (3.2) 
a 

I -  1)  ^  I  (Te  -  =  0  •  <5-3) 

Equations  (5.2)  and  (5.3)  are  identical  to  the  Hugoniot  conditions 
and  yield  the  following  discontinuous  solution: 


near  point  1  u^  =  u^  ,  ,  (5.^) 

near  point  2  u  =  u  *  1  ,  T  =  T»  =  1  .  (5.5) 

e.  d  ’  a  2 

This  discontinuous  solution  indicates  that  there  is  a  thinner 
atom  shock  imbedded  within  the  layer  of  thickness  t.  /c  and 
this  shock  will  be  studied  later.  Outside  the  shock  the  atom 
flow  is  uniform. 

From  Eqs.  (J».ll)  and  (l4,ll»)  we  find 


“l  '  “a  *  °''‘lay'll  > 


(5.6) 


T.  = 


(5.7) 


From  Fig.  2  we  recall  that  the  ratio  1.  /£..  is  very  small  and 

ia^ 

ct.  /I,,  is  even  much  smaller.  The  ion  slip  is  indeed  negligible, 
ia2  iig 
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Equatioo  (^.13)  yields  two  different  equations,  one  on  each 
side  of  the  discontinuity.  Near  the  point  b(b  *  1,2)  we  have 


-  dT 

3  e 

2 


e 


e 


u 

m 


-  T  )  {Q 
ra  ee 


(5.8) 


Equations  (ij.lO),  (1».9),  end  (l*.l6)  determine  respectively 
the  electric  ficla,  the  charge  separation  and  the  potential  by 
the  relations 


E 


6 


eXp 

2  dE 

t  3? 


* 


(5.9) 


(5.10) 


♦ 


(5.11) 


We  now  split  each  of  Eqs.  (5.8)  into  two  first  order  dif¬ 
ferential  equations  suitable  for  nximerical  integration.  The 
result  is 


t 


dT 

3  e 

2  dx^ 


.^1/2 


♦  b  -V  (T^  -  T  )  G 
<:  e  m 

u 

n 


(5.12) 


(5.13) 
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(5.11*) 


The  independent  variable  x*  is  next  eliminated  by  dividing  Eqs. 

(5.13)  by  Eq.  (5.12).  The  resulting  direction  field  equation 

=  f(z,  T^)  is  integrated  numerically  between  points  1  and  2 
e 

to  give  a  relation  z  =  z(T^).  With  z  known,  the  electron  tem¬ 
perature  distribution  T^(x*)  can  be  foxmd  by  integrating  Eq, 

(5.12)  where  the  origin  of  x*  is  taken  at  a  point  such  that 
-3 

+  10  .  The  electric  field,  charge  separation  and 

potential  distributions  are  computed  respectively  from  Eqs,  (5.9), 
(5.10),  and  (5.11). 

Ibe  upstream  and  downstream  points  1  and  2  are  obviously 
singular  points  of  the  system  of  Eqs.  (5.12)  and  (5.13)  and  the 
direction  field  equation  at  these  points  is  given  by  the  ratio 
of  two  vanishing  gradients.  It  is  therefore  necessary  to  determine 
the  nature  of  the  solution  in  the  neighborhood  of  the  singular 
points  by  a  suitable  linearization.  Since  at  the  singular  jjoints 
(m  =  1,2)  T^  =  T^  and  z  =  0,  we  assume  that  in  the  neighborhood  of 
the  singular  points  the  solution  is  of  the  form 


T  =  T  +  Ae 
e  m 


, 


(5.15) 


where  A  and  B  are  assumed  to  be  small.  On  inserting  these  assumed 
forms  of  the  solutions  in  Eqs,  (5.1?)  and  (5.13)  and  retaining 
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only  those  terms  linear  in  A  and  B,  one  obtains  a  homogeneous 
system  for  A  and  B.  The  condition  for  the  existence  of  a  non¬ 
trivial  solution  is  given  by  the  characteristic  equation.  The 
nature  of  the  roots  of  the  characteristic  equation,  here  a 
quadratic  equation  in  k,  determines  the  nature  of  the  singular 
points.  The  slope  of  the  characteristic  direction,  dz/dT^, 
at  these  points  is  simply  B/A. 

At  the  singular  points  m  =  1,2  the  ch«u*acteristic  equations 
of  the  system  of  Eqs,  (5.12)  and  (5.13)  are 


m 


=  0 


(5.16) 


where  F  and  G  denote  the  values  of  F  and  G  at  point  m. 
e  m  e 

m 

Equations  (5.16)  always  have  two  roots  of  opposite  sign.  It 


follows  that  points  1  and  2  are  saddle  points  and  the  characteristic 


directions  are  given  by  the  relation 


k  T^^^ 
dz  m  m 

dT  ”  F 

e  e 


(5.17) 


where  k^^  is  the  positive  root  and  k^  the  negative  one  of  the 
respective  characteristic  equations.  The  integral  curves  obtained 
by  integrating  the  direction  field  equation  from  points  1  and  2 
successively  are  sketched  in  Fig.  3. 
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It  has  been  pointed  out  that  there  is  a  disccxitinuity  in 

the  transition  from  point  1  to  point  2.  This  discontinuity  is 

represented  in  Fig,  3  by  the  line  3-**,  where  point  3  lies  on  the 

integral  curve  leaving  point  1  at  x*  ~  and  point  is  on  the 

integral  curve  reaching  point  2  at  x’-*-  +  ».  It  will  be  seen 

that  this  discontinuity  corresponds  to  an  atom  shock  of  thickness 

I  <<  £..  /e  in  which  the  ion  velocity  and  temperature  Jump 
8^2  ^^2 

from  their  upstream  to  their  downstream  values  but  where  the 
electron  temperature  can  be  considered  8is  constant.  In  particular 
we  have 


==  .  (5.18) 


Integrating  Eq.  (^4,13)  with  respect  to  x*  between  points 
and  3  with  the  assumption  of  constant  electron  temperature  we 
obtain  the  following  equation 


f(z3  -  zj^)  =  log^  u^  +  b 


(T  -  1)  G  dx«  . 


®3  ®3 


(5.19) 


The  function  G  depends  very  weakly  on  the  ion  velocity  and  may 
be  regarded  as  constant.  Since  the  shock  thickness  is  thin, 

X*  x|J  and  the  integral  on  the  right  hand  side  of  Eq.  (5.19) 
vanishes.  The  relation  between  the  coordinate  z  of  points  3 


and  becomes 
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z 


3 


(5.20) 


For  each  point  of  the  integral  curve  2-1*  ve  compute  the 
corresponding  point  3  from  Eq.  (5.20)  to  get  the  locus  of  point 
3,  Ihe  intersection  of  this  locus  and  of  the  integral  curve 
1-3  determines  the  actual  point  3.  The  solution  is  unique, 

Ihe  electron  temperature,  electric  field  and  potential 
distributions  in  the  thermal  layer  are  plotted  in  Fig.  h  and  5 

ll 

for  *  3  *  10  °K,  =  10  =  IQ-^  °K,  respectively 

-3  15  -3 

and  a  =  10  ,  n  =  10  cm  in  both  cases.  They  are  similar 
to  the  corresponding  distributions  obtained  in  the  fully  ionized 
case^.  Due  to  the  high  electron  thermal  conductivity,  the 
electron  temperature  rises  to  reach  a  value  close  to  its  down¬ 
stream  value  at  the  atom  shock.  Because  of  the  low  degree  of 
ionization  the  atom  gas  is  essentially  unaffected  by  the  electrons 
and  remains  in  its  upstream  state  ahead  of  the  shock  front.  In 
fact,  the  electrons  are  heated  up  by  a  very  small  decrease  of 
ion-atom  kinetic  energy  which  is  here  negligible  but  which  is 
appreciable  for  higher  values  of  a.  The  ion  flow  is  to  a  very 
good  approximation  in  equilibrium  with  the  atom  flow,  ht  =  3 
the  electron  temperature  rises  slowly  at  first  and  then  more  and 
more  rapidly  near  the  shock.  At  =  10  on  the  contrary,  the 
electron  temperature  rises  sharply  at  the  very  beginning  and  more 
slowly  afterwards.  After  the  shock,  the  electron  tem;p<*rature 
reaches  its  downstream  value  slowly  through  collisions  with  the 
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at>oss  and  ions*  The  thickness  of  the  theraaJL  layer  varies  froa 

31^^  /t  at  =  1.5  to  0.75  t..  /e  at  =  10. 

2  2 
The  electric  field  reflects  the  behavior  of  the  electron 

tenperature  distribution  since  it  is  proportional  to  its  gradient. 

It.  is  always  negative  and  for  high  l!ach  nuabers  exhibits  a  sharp 

single  oscillation  at  the  beginning  of  the  theraal  layer  as  a 

co;“Telary  to  the  sharp  rise  in  T^.  It  is  precisely  the  sane 

precursor  electric  shock  layer  which  has  been  discussed  ’n  the 

fuiUy  ionized  case^.  Its  thickness  is  M  I  .  /c  «  t.  /c.  The 

1  ii^  ii^ 

pot«‘ntial  distribution  is  identical  to  the  electron  temperature 

distribution  except  for  the  potential  Jump  across  the  shock. 

The  atom  shock  structure  equations  are  obtained  now  by 

sett:ng  A  =  £  ,  replacing  x^  by  x  =  x’/t  and  neglecting 

®  aa^  s  aa^ 

tenis  of  order  a  and  g.  Equations  (^*.13),  (I*. 12),  (*^.15),  (**.11), 
(*♦.13),  (J+.IO),  (*4.9),  and  (*4.l6)  yield  respectively  the  follow¬ 
ing  equations: 


dT  aa 

—  =0(ct^)  -0, 


(5.21) 


(hence  T  =  constant 
e 


t1/2 


du 


Q  dx 
aa 


T  ) 
3 


T 

1) .  -i 

a 


-  1 


(5.22) 


T^^^  dT 

^  Q  dx”  “  ^  ^  I  Q  dx  » 

aa  aa 


T^'^^  du 
a  a 


(5.23) 
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“i  *  “  a  i 


u.u 
i  a 


“a  *  ’'i> 
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T  +  T 

du.  _  _  i  e-  ,  dT. 

r— i.  (1  }?  _ h  +  1 _ 1 

^dx  S  S  “  2  ^  u.  dx 

u.  1 


^  du^ 

^  I  lx  ^*17"  dT"^  ^  • 
aa^  1 


(5.2U) 


T.  =  T  -  ^  MIu  (u.  -  u  )  -  ^-  ■■■■  . .  ■  - - -  v-f  ' 

1  a  9  2a  .  a  W  q  (j  * 

<;  ia  1  a 


^  dT^  T,  du^ 

^2  dx  *  ~  dx 
1 


^ia  du  dT. 

aa^  1  1 


(5.25) 


E 


du^ 

dx 


» 


D 

2  ^ 

t  dx  * 
aa„ 


log^ 


(^)  .  ^3 


(5.26) 


(5.27) 

(5.20) 


Equations  (5.22)  and  (5.23)  are  the  Navier-Stokes  equations  for  a 

neutral  gas  whose  solution  for  the  velocity  and  temperature 

17 

distributions  is  well  known  .  Since  1.  /£  <<  1  Eqs.  (5.2^4) 

ia_  aa„ 


17.  D,  Gilbarg  and  D.  Paolucci,  J.  Rat.  Mech.  and  Anal,  2,  6l7, 
(1953). 
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and  (5.25)  show  that  to  a  good  approxination ,  u.  =  u  emd 

X  & 

T.  =  T  except  at  the  beginning  of  strong  shocks  where  the 

X  & 

dimensionless  temperatures  are  very  small.  Equation  (5.28) 
shows  that  the  electrons  in  thermal  equilibrium  obey  a 
Boltzmann  distribution 


n'  =  n'  exp  (e(({(‘  -  j 

3  ^  3 

The  atom  velocity  and  temperature  distributions  are  obtained 
by' solving  the  system  of  Eqs,  (5.22)  and  (5.23)  according  to  the 
procedure  used  earlier  to'  integrate  the  electron  temperature 
distribution.  Points  3  and  are  singular  points  of  this  system 
and  the  characteristic  equation  at  point  h  is: 

dck^  +  k(d(l  -  I  J^)  -  I  c)  +  I  (M^  -  1)  =  0  .  (5.29) 


Since  <  1,  Eq.  (5.29)  has  two  roots  of  opposite  sign.  Point 
is  a  saddle  point.  The  characteristic  direction  of  the  integral 
curve  in  the  T  -  u  plane  is  given  by 

Gl  GL 


dT 

- 

a 


(5.30) 


tfhere  k  is  the  negative  root  of  Eq.  (5.29). 

The  characteristic  equation  at  point  3  is  deduced  from  Eq, 
(5.29)  by  replacing  by  in  the  parameters.  Since  >  1,  it 
is  easy  to  verifj'’  that  the  characteristic  equation  at  point  3  has 
2  positive  roots  and  point  3  is  a  node. 
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Ijivrre  is  only  one  integral  curve  passing  through  the  saddle 
point  with  the  characteristic  direction  given  by  Eq,  (5.30), 
while  the  node  3  is  a  center  of  attraction  for  edJ.  integral 
curves  passing  in  its  neighborhood.  Therefore,  the  numerical  in¬ 
tegration  is  performed  from  point  ^  towards  point  3  and  the  origin 
of  X  is  chosen  at  an  arbitrary  point  near  point  U. 

When  the  ion  slip  is  small  the  ion  velocity  and  temperature 

distributions  are  assumed  to  be  given  by  the  following  expansions 

in  the  small  parameter  t.  ft 

ia^  aa^ 


(5.31) 


The  terms  u.  , ♦ , ,  T.  ...  are  computed  by  substituting  the 

h  ^1 

expansions  (5.31)  in  Eqs,  (5.2*4)  and  (5.25)  and  identifying  the 

terms  of  t  .e  same  order  in  £ ,  /£  .  The  first  order  terms  are 

ia^  asg 


u, 


2u^  dn  dT 

- - - T7?  - P - ^5.32) 

aoTurT^  dx  3  2  ^2  u^d. 

ia  a  a 


*4U‘= 


dT  T  du 


[■  »-|m?uu - 2 - .  (5.33) 

^1  9  2  a  i^  ^  ^ 

X  &  8t 


We  note  that  the  first  order  approximation  for  u^  and  T^  does  not 
contain  the  effect  of  ion  viscosity  and  thermal  conduction  which 


are  very  small 


The  electric  field,  potential,  velocity,  and  temperature 
distributions  for  atoms  and  ions  are  plotted  in  Figs,  6  and  7 
for  =  3  and  10  respectively.  The  atom  flow,  unaffected  by  the 
charged  particles,  undergoes  a  shock  in  which  the  dissipative 
mechanisms  are  the  atom  viscosity  and  thermal  conduction.  Be- 

S 

cause  of  the  large  atom-ion  cros;8^section,  the  ion  flow  follows 
closely  the  atom  flow  and  goes  through  the  same  shock.  The 
negative  electric  field  induced  by  the  charge  separation  slows 
down  the  ions  so  that  the  ion  velocity,  which  is  at  first  slightly 
larger  than  the  atom  velocity,  then  decreases  more  rapidly.  The 
ion  temperature  increases  faster  than  the  atom  temperature  and 
the  atom-ion  equilibrium  is  reached  at  the  end  of  the  shock  after 
a  sufficient  number  of  atom-ion  collisions. 


At  =  3  the  ion  slip  is  small  throughout  the  shock.  At 
=  10,  the  ion  slip  is  still  small  except  at  the  beginning  of  the 
shock  where  the  method  of  solution  based  on  the  expansions  (5.31) 


fails,  as  may  be  seen  on  Fig.  7.  No  acceptable  numerical  solution 

for  u^  and  T^  in  the  neighborhood  of  point  3  has  been  found  yet. 

1/2 

Analytically,  the  solution  fails  because  T  is  very  small  near 

Bi 

point  3  and  the  first  order  terms  given  by  Eqs,  (5.32)  and  (5.33) 

are  lairge.  A  physical  explanation  might  be  that  the  local  collision 

—1  /2 

time  proportional  to  T"  '  is  much  larger  at  the  beginning  of  the 


shock  than  at  the  end  because  of  the  high  temperature  ratio  across 
the  shock.  As  a  result  there  are  not  enough  ion-atom  collisions  at 


-loa¬ 


the  beginning  of  the  shock  to  keep  the  ion  slip  small.  A  smaller 
value  of  the  ratio  I,  /l  would  improve  the  range  of  validity 

lag  aag 

of  our  solution. 

The  electric  field,  always  negative,  shows  a  single 
symmetrical  oscillation.  The  ratio  R  of  the  Coulomb  force  acting 
on  the  ions  to  the  ion  pressure  gradient,  defined  by  the  relation 


dx* 


increases  through  the  shock  from  Q^k  to  0.8  for  »  10.  There¬ 
fore,  there  is  an  Important  coupling  of  the  electrical  effects 
with  the  ion  flow. 

The  magnitudes  of  the  shock  thicknesses  and  of  the  electrical 
effects  are  given  in  Table  1  for  typical  plasma  conditions  at 
different  Mach  numbers.  The  thickness  of  the  thermal  layer  ahead 
of  the  shock  front  is  denoted  by  The  shock  thickness  based 

on  the  maximum  atom  velocity  gradient  is  denoted  by  Ag  and  is  given 
by 


^1  ■  ^ 
du 

(—) 

dx  meue 


As  the  Mach  number  increases  the  dimensionless  electric  field  and 
the  potential  increase  while  the  shock  thickness  increases.  For 
different  temperature  and  density  conditions  the  electric  field  E* 

3/14 

would  vai*y  as  n'T’  and  the  potential  rise  as  Ti, 
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VI.  PARTIALLY  IONIZED  PLASMA 
We  now  consider  in  detail  the  case  in  which  the  degree  of 
ionization  a  is  no  longer  small.  In  this  case  the  ionized 
particles  have  an  important  effect  on  the  atom  flow  and  the 
magnitude  of  this  effect  increases  with  the  degree  of  ionization. 
We  wish  to  investigate  how  the  shock  structure  and  the  self- 
induced  electric  field  vary  with  the  upstream  conditions  and  the 
degree  of  ionization.  We  shall  compare  the  results  obtained  with 
the  limiting  case  of  a  fully  ionized  plasma. 

The  shock  structure  is  now  investigated  according  to  the 
method  employed  previously.  We  consider  fxrst  the  thermal  layer 
of  thickness  A  =  £..  /e  where  x  is  replaced  by  x*  =  x/l..  . 

s  lij  S  iig 

Equations  (5,6)  and  (5,7)  still  hold  and  we  may  set  u^  =  u^, 

T.  =  T  in  the  remaining  shock  equations.  On  neglecting  those 
X  s, 

terms  of  order  e  but  retaining  those  of  order  a  in  Eqs,  (^^,12), 
(^.15)t  (^,13),  (^,10),  (^,9),  and  (^,l6)  respectively  we  obtain 
the  following  shock  structure  equations: 

T  T 

I  (1  ♦  a)  (u^  -  1)  +  ~  -  1  +  a  (^  -  1)  -  0  ,  (6,1) 

a  a 


Parr'^^  dT 

e  e 

5F  dx^ 
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(l  +  a)(u‘'-l)  +  T  -  l-^a(T  -  l) 
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(6,2) 
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where  G  is  defined  by  (5.1^), 


p  2  d  f  e. 

hi^  “a 


eX, 
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ii2  r 
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E  dx» 


f 

=  u  d(-^) 
/  a  u 

^  ft 


(6.i») 


(6.5) 


(6.6) 


Equation  (6,1)  nay  be  solved  for  T  and  the  resulting  relation 

& 

used  to  eliminate  T  in  Eq.  (6,2)  to  give 

A 


dT^  lOM^d  +  c) 

^  "F“  dx’f  “  3^ 

e 


(6.7) 


where  u^^  *  (3  +  . 

Differentiating  both  sides  of  Eqs,  (6,7)  with  respect  to  x*, 

substracting  the  resulting  equation  from  Eq.  (6,3)  and  solving 
du 

for  r— y  we  obtain 
dx* 


du 

a 

d?r 


-  T^)G][ 


^(1 


2a  ^  3  2  a  3  2 


+  1) 


(6.8) 


Once  T  is  eliminated  in  Eq.  (6,8)  with  the  aid  of  Eq,  (6.1), 

Gl 

Eqs,  (6.7)  and  (6,8)  form  a  set  of  two  first  order  differential 


equations  in  T  and  u  .  The  method  of  solution  is  the  same  as 

6  S' 

that  described  earlier  for  the  veakly  ionized  case.  The  variable 
X*  is  eliminated  by  dividing  Eq.  (6,8)  by  Eq.  (6.7)  and  the 
direction  field  equation  du  /dT  =  f(T  ,  u  )  is  integrated 

S  S  6  S 

between  points  1  and  2  to  give  a  relation  u^  =  '^a^'^e^* 
electron  temperature  distribution  T^(x*)  is  obtained  by  integrating 
Eq.  (6,7)  and  the  ion  temperature  distribution,  electric  field, 
charge  separation,  potential  are  given  in  turn  by  Eqs.  (6.1), 

(6.^4),  (6.5)  end  (6,6)  respectively. 

Points  1  and  2  are  again  singular  points  of  Eqs,  (6.7)  and 
(6.8)  and  the  characteristic  equation  at  the  points  m  =  1,2  is 


fk^  .  2a 

F  S(l  +  a) 
e 
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1  +  wf)  k(|  (1 
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m!)  .  -  1)J 
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m 


+  b  G  (1  -  M^)  (1  +  a)  =  0 
m  m 


(6.9) 


Ihe  characteristic  directions  are 


2afk 

5(1  +  a)F  (1  ~  M^) 
e  m 


(6.10) 


m 

The  product  of  the  roots  of  Eq.  (6.9)  is  of  the  sign  of  the 

expression  (1  -  M^)  ( ^  point  1, 

m  5v  -L  +  a ;  m 

O  Of* 

f-C  >  1  >  1  -  r/vv”  T  the  product  of  the  roots  is  negative, 
1  5v  1  Q/ 

Eq.  (6.9)  has  2  roots  of  opposite  sign  and  point  1  is  a  saddle 
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point.  At  point  2  tvo  different  situations  arise:  if 

>  1  -  Eq.  (6.9)  has  two  negative  roots  and  point 

2  is  a  node,  The  mininum  value  of  is  reached  at  a  =  1 
(fully  ionized  plasma)  and  corresponds  to  a  maximum  free  stream 


Mach  number  of  1,12.  This  situation,  therefore,  corresponds 
to  weak  shocks.  If  <  1  -  ,  Eq.  (6,9)  has  two 

roots  of  opposite  sign  and  point  2  is  a  saddle  point. 

As  in  the  fully  ionized  case^  but  in  contrast  to  the 
weakly  ionized  case,  fcr  a  #  0  there  are  two  different  types 


of  shock  structure  depending  upon  the  Mach  number.  The 
particular  "weak  shock"  solution  cannot  exist  in  the  weakly 


ionized  case  because  the  condition  1  - 


be  satisfied  when  a  ^  0, 

2 

In  the  weak  shock  case,  when  M^  >  1 


r  <  M?  <  1  cannot 
•♦■a)  2 


,  the 


numerical  integration  is  carried  out  from  point  1  to  point  2. 


The  solution  in  the  u  ,  T  phase  space  is  continuous  and  has  no 

extrema.  The  shock  structure  consists  of  a  single  relaxation 

layer.  The  results  for  =  l,07l<  at  a  =  0,7  are  plotted  in 

Fig,  8.  The  shock  thickness  is  so  lEirge  (£.  /c)  that  the  atoms, 

li2 

ions,  and  electrons  make  many  collisions  and  are  in  mechanical 
and  thermal  equilibrium  through  the  shock. 

When  ^^  <  1  -  integral  curve  pattern  in  the 

u^,  Tg  phase  space  is  sketched  in  Fig.  9.  The  integral  curve  1-3 
leaving  point  1  with  the  appropriate  characteristic  direction 
(corresponding  to  a  positive  root  of  Eq,  (6.9))  has  a  maximum  at 


point  5  where  it  intersects  the  isocline  of  zero  slope,  a  parabola 
whose  equation  is 


T 

e 


+  aj 


,  u  (.  I  .  I 


3  2  a  3  '2 


1) 


♦ 


and  does  not  reach  point  2  with  the  proper  characteristic  direction 
(corresponding  to  a  negative  root  of  Eq.  (6.9)).  Similarly,  the 
integral  curve  2^h  leaving  point  2  with  the  appropriate  direction 
does  not  reach  point  1.  There  is  no  continuous  integral  curve 
Joining  points  1  and  2  and  we  must  adopt  the  discontinuous  solution 
l_3_b-2.  The  discontinuity  between  points  3  and  h  represents  a 
shock  layer  much  thinner  than  the  thermal  layer  which  is  investigated 
below. 

V/e  know  that  the  shock  thickness  is  of  the  order  of  Z  for 


aa^ 

a  weakly  ionized  plasma  and  of  the  order  of  t..  for  a  fully  ionized 

ii2 

plasma.  We  may  expect  that,  for  intermediate  values  of  the  degree 
of  ionization,  the  shock  thickness  is  a  function  of  the  three  mean 


free  paths  £  .  Considering  the  case  in  which  a  is 

an2  la^ 

small  (j  0.3)  we  again  set  A  =  £  .  Eq.  (^,15)  gives: 

s  aa^ 


dT 

e 

dx 


0(  c 


£ 


)  .0 


(6.11) 


Hence,  the  electron  temperature  T^  may  be  considered  as  constant 
inside  the  shock  and  we  set 


T 

e 


_  Y 


(6.12) 


-Ii8- 


On  using  Eq.  (6.12),  Eqs,  and  {h,lU)  yield 


u.  «  u  - 

i  a 


2  t  u.u 
ia^  i  a 


du. 


T  +  T 
i  e. 


at  Q,  (T  T. ) 

aa2  ia  a  i 


m  'dT  <1 4^^  *  “>  -  2 


u. 


i-fi. 

u,  dx“^l  dx  F.  dx^^  • 
i  aa^  i 


(6.13) 


Ul  u  u 

c  p  2.ap  i  a 

■^i  •  '^a  “  I  ^  V'^i  “  - ^ - 172 

^  a  9  ^  a  i  a  bl  Q,  (T.  +  T 

aa^  ia  i  a 


^2  dx  u,  dx  “  1  F,  ^dx  ^ 

i  aa^  i 


dT 

*^lr<4;-dr”  • 


(6.14) 


The  ion  slip  is  of  the  order  of  1.  1 1  ,  therefore  small.  It 

la^  aa2 

follovs  that,  for  a  small,  we  may,  to  a  good  approximation  replace 

u,  by  u  and  T,  by  T  in  Eq.  (^*.12)  which  becomes 
la  la 


^t1/2  %  (1^  ,  IS,  .  5  ^,1  ,  „)(„  .  1,  ,  !a 

dx  F  F.  I  3  2  a  u 

a  i  aa^  a 


a 


T 

1  +  a(^  -  1) 


(6.15) 


It  is  convenient  at  this  point  to  replace  Eq.  (^4.13)  by  the 
equation  obtained  by  adding  the  atom  and  ion  energy  equations 
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(2.6)  which,  in  dimensionless  form,  is  written 


a  1  a* 


(1  +  o)(l  +  f  '4>  “a  *  I  ’'a  ■  “a*  '  ’ 


(6.16) 


The  derivation  of  Eq.  (6.l6)  is  given  in  the  Appendix  B,  The 
constant  of  integration  F  is  determined  at  the  boundary  conditions 
at  point  3  or  1*. 

The  Hugoniot  relations  for  the  shock  provide  the  boundary 
conditions  at  points  3  and  h  and  are  obtained  by  setting  the 
derivatives  in  Eqs,  (6,15)  and  (6.l6)  equal  to  zero.  They  are 


I  1^(1*  a) 


-l)+-2--l-faT  =0  , 

u  e_  * 

a  3 


+  4  T  -  aT  log  u  -  F  =  0 
2  a  e^  e  a 


(6,17) 


(6.18) 


The  elimination  of  T  between  Eqs,  (6,17)  and  (6,l8)  and  of  F 

a 

computed  at  point  U  gives  the  "Jump”  condition  -satisfied  by 
the  velocity  at  point  3 


5(|  +  l)(u^  r 


u 


c  u, 


=  0 


(6.19) 
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If  ve  put 


y 


1 


2a 

1  +  a 


u 


and 


then  from  Fig.  10  it  is  clear  that  the  equation  equivalent 

to  Eq,  (6.19),  has  only  one  root  u  other  than  u  *  u,  ,  For  each 

point  of  the  integral  curve  2-k,  we  compute  the  corresponding  jxsint 

3  from  the  Jximp  condition  (6.19)  to  get  the  locus  of  point  3.  The 

intersection  of  this  locus  and  of  the  integral  curve  1-3  leaving 

point  1  determines  the  actual  point  3.  The  solution  is  unique. 

Ibe  structure  equations  (6.15)  and  (6,l6)  are  integrated  in 

the  same  way  as  the  thermal  layer  equations.  The  characteristic 

equation  at  the  singular  points  m  »  3,^*  obtained  by  setting 
lex 

u  «  u  ♦  Ae  ,  etc. ,  is  here 


T  (1  -  a)  T^''^k  ’^m  “'^e 

dc(l  -  a)^  J - - 2 -  [d(|  iqa  ♦  a) - 5 - ^) 


u 


m 


( 3T  ♦  «T  ) 

♦  J^(l  ♦  a)  -  (|  M?  ♦  1) - - - ^  .  0  ,  (6.20) 

M  2u‘^ 


where  F  is  the  value  of  the  function  F  at  point  m.  Point  3 

ft  ft 

m 

turns  out  to  be  a  node  aud  point  a  saddle  point.  The  atom 
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velocity  and  temperature  distributions  are  obtained  by  integrating 
the  direction  field  equation  dT  /du  =  g(u  ,  T  )  from  point  to 
point  3. 

The  ion  velocity  and  temperature  distributions  may  not  be 

given  by  the  expansions  (5.31)  when  the  shock  thickness  is  smaller 

than  1  .  However,  retaining  the  assvmption  of  small  ion  slip 

aa^ 

we  substitute  u  and  T  for  u.  and  T,  in  the  right  hand  side  of 
a  a  1  1 

Eqs.  (6.13)  and  (6.ll*).  This  iterative  method  yields  the  follow¬ 
ing  first  order  approximations  for  u^  and  T^ : 


u. 


a  - 

a 


21,  u" 

ia^  a 


al  Q.  (2T  ) 
aa^  la  a 


i72  '3r  ‘f  - 


T  +  T 
a  e- 


u 


•f 


u  dx 


a 


du 


(6.21) 


^1  *  \  “a<“i  -  V  - 


l4t.  U 

ia„  a 


bt  Q,  (2T  ) 
aa^  ia  a 


dT 

2  r3  a 

172  2  dx 


T  du 

a  a 

u  dx 
a 


S 


<— >  *  ST ' 


t1/2 


a 

F, 


dT 

a 

dx 


■)] 


(6.22) 


It  is  found  that  the  ion  viscosity  and  thermal  conduction  appreciably 
reduce  the  ion  slip  and  must  be  taken  into  account  for  a  >  0.1. 
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As  in  the  weakly  ionized  case,  the  electric  field,  charge 
separation  and  potentiail  are  given  respectively  by  Eqs.  (5.26), 
(5.27).  and  (5.28). 

Results  for  two  different  degrees  of  ionization  a  =  0.1 
and  a  =  0.5  at  a  free  stream  Mach  number  =  2  are  presented 
in  Figs.  11  and  12  respectively.  In  contrast  to  the  weakly 
ionized  case,  the  atom  shock  retains  some  features  of  the  icn 
shock  in  the  fully  ionized  case^. 

In  the  thermal  layer  the  electrons  are  heated  up  by  the 
decrease  in  the  kinetic  energy  of  the  flow,  faster  than  the  heavy 
particles,  because  of  their  higher  thermeuL  conductivity.  When 
the  degree  of  ionization  becomes  higher^  the  thermal  energy  of 
the  electrons  represents  a  larger  fraction  of  the  total  energy 
and  the  velocity  drop  and  the  heating  of  atoms  and  ions  increase. 

The  heating  and  the  compression  of  the  heavy  particles  occur, 
mainly  in  the  shock  layer  where  it  can  be  seen  that  the  atom  and 
the  ion  temperatures  overshoot  their  downstream  values.  Because 
of  their  much  smaller  mass  the  electrons  do  not  have  time  in  a 
few  mean  free  paths  to  exchange  energy  with  the  heavy  particles 
and  the  electron  temperature  is  approximately  constant  through 
the  shock. 

The  dissipative  mechanisms  in  the  shock  are  the  viscosity 
and  the  thermal  conduction  of  the  ions  and  atoms.  The  ion  slip 
is  everywhere  small  for  a  =  0.1;  it  is  considerably  larger  for 
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a  =  0.5.  Moreover  it  may  be  seen  on  Fip.  12  that  the  ion  tem¬ 
perature  and  velocity  distributions  obtained  from  Eqs.  (6,21) 
and  (6.22^  are  locally  in-  error  at  both  ends  of  the  shock  where 
they  overshoot  their  boundary  values.  In  fact,  for  a  =  0.5  the 

solutions  u  (x)  and  T  (x)  obtained  from  Eqs.  (6.15)  and  (6.I6) 
a  a 

do  not  represent  the  true  atom  velocity  and  temperature 
distributions  but  rather  the  average  distributions  for  the  whole 
plsisma.  Ihis  approximation  may  be  responsible  for  the  erroneous 
behavior  of  the  ion  velocity  and  temperature  distributions  at 
both  ends  of  the  shock.  Unfortunately  no  acceptable  numerical 
solution  which  permits  a  simultaneous  calculation  of  the  ion  and 
atom  velocity  and  temperature  distributions  has  as  yet  been  found. 

The  electric  field  varies  slowly  in  the  thermal  layer  and 
there  is  a  single  oscillation  of  negative  amplitude  in  the  shock. 
The  potential  reaches  its  maximum  value  at  downstream  infinity. 

The  shock  structure  for  a  strong  shock  (M^  =  10)  at  a  =  0.1 
is  shown  in  Fig.  13,  As  in  the  weakly  ionised  case  the  electron 
temperature  rises  sharply  at  first  and  the  electric  field  has  a 
very  sharp  negative  oscillation  with  a  slower  damping  in  the 
precursor  electric  shock  layer  at  the  beginning  of  the  thermal 
layer.  The  ion  slip  is  more  important  than  at  »  2  for  the  same 
degree  of  ionization  and  for  the  same  reason  as  in  the  weakly 
ionized  case  the  ion  velocity  and  temperature  distributions 
obtained  from  Eqs.  (6.21)  and  (6.22)  are  in  error  near  point  3. 

(The  erroneous  parts  are  drawn  in  dotted  lines  on  Figs.  12  and  13.) 


The  variation  of  the  shock  thickness  A  and  of  the  jsaxiaaa 

s 

electric  field  and  potential  with  the  degree  of  ionization  a  is 
given  in  Table  II  for  =  2  and  =  10,  When  a  increases,  the 
total  number  of  atoms  and  iems  being  constant,  the  maximum  electric 
field  measured  in  V/cm  increases,  but  since  the  shock  thickness 
{meas^lred  in  cm)  decreases  the  potential  rise  aci^css  the  shock 
stays  constant  for  all  degrees  of  ionization.  The  shock  thickness 
is  a  minimum  and  the  electric  field  is  a  maximum  when  the  plasma 
is  fully  ionized.  The  shock  thickness  is  smaller  than  the  down¬ 
stream  atom-atom  mean  free  path  t  for  o  >  0.2  but  is  always 

aa2 

much  larger  than  the  downstream  ion-ion  mean  free  path  I  . 

ii2 

except  when  a  =  1,  The  shock  thickness  decreases  with  the  Kach 
number  while  the  electric  field  and  potential  increase, 

VII.  CONCLUDIHG  REHARKS  AND  COMPARISON  WITH  EXPERIMENT 
The  model  that  we  have  proposed  covers  in  principle  all  plasma 
regimes  from  weakly  to  fully  ionized.  But  unfortunately,  the 
difficulty  of  numerically  integrating  more  than  two  simultaneous 
differential  equations  between  singular  points  requires  that  the 
solution  be  limited  to  an  ionization  less  than  about  30%  so  that, 
the  ion  slip  s  's  small  throughout  the  shock. 

Some  feacures  of  the  shock  structure  are  common  to  all  plasmas 
regimes  such  as  the  broad  region  of  elevated  electron  temperature 
ahead  of  the  shock  fre-nt  and  the  induced  electric  field  which  travels 
with  the  shock  front.  The  potential  rise  across  the  shock  is  in¬ 
dependent  of  the  degree  of  ionization. 
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When  the  degree  of  ionization  a  varies  froia  very  low  values 
to  one,  the  shock  strength  decreases,  and  a  partial  compression 
and  heating  of  the  heavy  particles  occurs  ahead  of  the  shock 
front.  The  atom  and  ion  temperatures  overshoot  their  downstream 
values.  The  shock  thickness,  which  is  of  the  order  of  the  down¬ 
stream  atom-atom  mean  free  path  when  a  «  1  and  of  the  order  of 
the  downstream  ion-ion  mean  free  path  when  a  «  0(l),  decreases 
continuously  when  a  goes  from  zero  to  one. 

The  potential  rise  across  the  imbedded  shock  in  a  magnetic 
annular  shock  tube  was  measured  by  Heywood  in  partially  ionized 
hydrogen  at  a  free  stream  temperature  of  300°K,  a  molecular 
density  of  8,75  •  10  cm  and  a  shock  speed  of  12,5  cm/ys. 

An  electric  probe,  consisting  of  two  steel  electrodes  separated 
by  a  distance  of  0,6  cm  and  insulated  from  each  other,  was  placed 
parallel  to  the  axis  of  the  shock  tube,  half  way  between  the  in¬ 
ternal  and  external  diameters,  Ibe  electrodes  were  connected  to 
the  poles  of  an  oscilloscope  which  measured  the  axial  potential 
difference  between  them  (voltage).  The  voltage  distribution  with 
time  is  shown  on  Fig.  lU.  As  the  travelling  shock  reaches  the  first 
electrode,  the  voltage,  initially  zero,  increases  to  a  sharp  peak  of 
36V  and  decreases  to  zero  again.  The  voltage  variation  occurs  over 
a  travel  tLme  corresponding  to  a  distance  of  1  cm.  (ihe  voltage  in¬ 
crease  on  the  right  is  due  to  the  driving  currents  of  magnetic 
origin  and  is  irrelevant  here.)  This  result  shows  that  the  shock 


l8.  J.  B.  Heywood,  Ph.D.  Thesis,  Department  of  Mechanical  Engineer¬ 
ing,  M.I.T.,  Cambridge,  September  196i».  To  be  published. 
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thickness  is  at  nost  equal  to  1  cn  al  .hough  not  smaller  than  the 
distance  between  the  electrodes}  were  it  otherwise,  the  voltage 
curve  would  have  a  flat  maximusn.  Therefore  in  this  case,  the 
shock  thickness  is  of  the  order  of  0.6  cm,  the  distance  between 
the  electrodes.  It  may  be  possible  in  future  experiments  to 
evaluate  the  shock  thickness  more  accurately  by  varying  the 
distance  between  the  electrodes.  The  maximum  voltage  occurs  when 
the  shock  is  exactly  between  the  electrodes  and  represents  the 
potential  rise  across  the  imbedded  shock.  For  the  conditions  of 
the  experiment,  assuming  that  the  dissociation  is  complete  and 
of  the  hydrogen  is  ionized,  which  corresponds  to  a  free  stream 
Mach  number  of  60,  the  theory  predicts  a  shock  thickness  of  0.25  can 
and  a  potential  rise  across  the  imbedded  shock  of  33V.  The  agree¬ 
ment  with  the  messiured  potential  rise  of  36V  is  excellent,  the 
agreement  with  the  approximate  measurement  of  the  shock  thickness 
is  also  good.  These  results  suggest  that  the  measurement  of  the 
potential  rise  may  actually  be  a  relatively  simple  method  for  the 
determination  of  shock  thickness  and  the  shock  front  location  in 


shock  tube  experiments 
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APPEHDIX  A 

ATOM-ION  ENERGY  TRMSFER 

Following  Byron  and  Petschek^^  the  ion -at cm  energy  transfer 
per  unit  volume  due  to  the  random  motion  of  the  particles  is 

HI  •  f 

AE  =  —  n’n^  1(1*^)  cj(g)  g  f^f^  *  (^.1) 

where  g  is  the  relative  velocity  in  the  mass  center  system, 

G  the  center  of  mass  velocity,  o{g)  the  cross  section  for  ion- 

atom  collisions,  f  and  f.  are  arbitrary  normalized  distribution 

a  1 

functions  and  c.  and  c  the  random  particle  velocities.  The  atom 
la 

mass  is  taken  equal  to  the  ion  mass  m^. 

We  choose  c  and  g  as  variables  of  integration.  According  to 

A 

the  equations 


Ci  «  g  Ca  . 


(A. 2) 


dc^  »  dg  , 


(A. 3) 


c  c. 

-►  /  a  i\ 

g  .  G  *  g  .  ( - - - ) 


c  •  g  s*”  . 

a  2  * 


(A.l*) 


Eq,  (A.l)  becomes 


AE 


2^  "K  I  «  Vi  ‘‘‘a  *  IT  I  Vi 


(A. 5) 
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Substituting  into  Eq.  (A, 5)  the  Maxwell  distributions 


m. 


•  a 


(A.6) 


^iCCa  +  g)  =  (2TfkT!^  21^  ^ 

1  1 


we  obtain 


m. 

AE  =  7^  n!n'  (• 
2  1  a 


^i  ^  r  ^i^  f 

; - ; - r  {  o(g)  g  exp  (- 

2Trk  /fTf^  J  i  J 

1  a 


m.c  , 

-*■  ■*•  r  1  a  /I 

=a  •  8  “P 

1 


+  —  >•  Va 
a 


g/kT’]  dc^  dj  +  i 


g  o(g)  exp  (-  r^) 


2kT! 

1 


•  f  axp  1-  ^  (iy  ♦ 


1  a 


m.c  •  g/kTI ]  dc  dg)  (A. 8) 

1  a  Of  a 


Integrating  Eq.  (A, 8)  with  respect  to  the  atom  velocity  and 
transforming  the  integration  over  the  velocity  g  into  an  integration 
over  the  speed  g  we  obtain 


®'i  2/2  *^a  5  **^i^ 

AE  =  27rmj  n^n^  ^2iTk(T’  +  T’ )  ^  T!  v  T*  ^  ^T!  +  T' 

a  i  1  a  ^  i  a 

«  2 
if  5 

^  g  o(g)  exp  [-  2'k(T*"  ■>•  t'M^ 

i  i  a 


dg 


(2k) 


1  ^/P  f  K 

.  x3/2  "a*^!  ^T»  V  T!^  “  "^a^  j  ®  2k(T!  +  T'T^  * 

JC  /  &  X  1ft 


(A. 9) 


Assuming  that  the  cross  section  a(g)  is  independent  of  g  as  in  the 
hard  sphere  model  the  integral  term  is  easily  computed  and 
the  result  is 

AE  *  (T*  +  Q?  njn*  k(T!  -  T» )  ,  (A.IO) 

/-  m.  i  a  ia  i  a  i  a  * 

i 

where  ve  have  set  o(g)  =  Q!  , 

la 

Hie  total  energy  transfer  is  the  simi  of  the  energy  transfers 
due  to  the  random  and  directed  motions  of  the  particles  and  is 
written 
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4. -------- HI, (u’  +  uO]  .  (A. 11) 
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APPENDIX  B 

ION~ATOM  ENERGY  EQUATION 

On  adding  the  ion  energy  equation  (2,6,  J  ■  l)  to  the  atom 
energy  equation  (2.6,  j  =  a),  we  obtain  the  following  equation 

C  du'^  du!^  c  , 

0“  j""  0“  “TT  o  Tnr  (n'u*  kT’  n'u!  kTj )  -  eE* 

1  2  dx*  i  2  dx'  2  dx'  a  a  a  i  i  i  i  x 


d 


du'  ,  du’  ,  dT'  ,  dT‘ 

(u'u”  “  TT  TT)  “  TT  Tt)  “  TT  TT^ 

a  a  dx'  dx'  i  i  dx'  dx’  a  dx'  dx'  i  dx' 


=  ^  +  u'P  +  +  u!P. 

ae  a  ae  ie  i  ie 


(B.l) 
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Inside  the  shock  of  thickness  I  or  ,  the  momentum  and 
energy  transfer  between  the  electrons  and  the  heavy  particles  is 
negligible  and  the  right  hand  side  of  Eq.  (B.l)  vanishes.  Further¬ 
more,  inside  the  shock  is  constant  and  Eq.,  (2,10)  yields 


C.  kT*  du! 

C.  eE’  =  "T" 
IX  u!  dx' 


(B,2) 


On  substituting  C^eE^  by  its  value  from  Eq,  (B,2)  in  Eq, 
IJ,1)  and  integrating  with  respect  to  x’  when  is  constant, 
we  obtain  the  following  equation 


dT*  du' 
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K  r  *  u'p  — r 
a  dx'  a  a  dx' 
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dT» 
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du* 

+  u.uV 
X  I  dx' 


2  S 

+  m,  ~  u!'"  ♦  k(C  T'  +  C.T')  -  C.  kT'  log^  u'  -  F'  ,  (B,3) 

1  d  1  d  aa  11  le  ei 


which  in  dimensionless  form  may  be  written 
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In  the  shock  we  set  A  »  t  ,  x  =  x,  u  »  u^,  T  =  T  ,  T  =  T 

8  ®  ^  a’  i  a*  e  e^ 

and  Eq,  (B,M  becomes 
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al. 


a  i  aa2 


*  I  "^a  -  “  ’’e  I”®®  “a  -  '■  • 


(B.5) 


Multiplying  Eq.  (6.15)  by  and  subtracting  from  the  above 
equation,  ve  obtain  Eq.  (6.I6). 


TABLE  2.  Maximum  maf^nltudes  of  electric  field,  potential  rise  acroaa  thock  and  ihock  thicknesi 
function  of  the  degree  of  ionization  a  for  moderate  and  strong  shocks. 
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POTENTIAL  <f> 


ELECTRIC  FIELD  E 


Electron  temperature,  electric  field  and  potential 
distributions  In  the  thermal  layer  for  the  weakly 
Ionized  case  (c  *  lo"^)  at  M.  =3. 
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FiR.  5.  Electron  temperature,  electric  field  and  potential  distributions  in  the  thermal 


layer  for  the  weakly  ionized  case  (o  «  10“^)  at  =  10. 
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Velocity,  temperature,  electric  field  and  potential 
jiistributions  in  the  atom  shocH  for  the  weakly 

—3 

ionized  case  (a  =  10  )  at  =  3* 
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Velocity,  temperature,  electric  field  and  potential  distributions  throu^ 
a  weak  shock  (M^  ■  1.07^)  at  a  «  0.7. 
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Fig.  13.  Velocity,  temperature,  electric  field  and  potential  distributions 
through  a  strong  shock  =  10)  at  a  =  0,1, 
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Fig.  Voltp^  distrlb^ltioc  as  a  ^f^ction  of  tine  for  a  shock 

spe^  of  12.5  or/vs,  %  free  str«^  tai^rature  of  300^K 
and  a  siolecular  densltj  of  8,75  •  10^^  c*"^  in  partially 
iooixed  hydrogen. 


Fig.  ll*b. 


Detail  of  the  region 


shown  by  the  arrow  in  Fig.  ll*a 


corresponding  to  the  potaitial  rise  across  the  inbedded  shock 


